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Abstract 



> 

►7- . Some atomic decays of elementary particles with pionium (7r + 7r~-dimesoatom) in 

the final state are considered. Only for K-meson atomic decays are the corresponding 
branching ratios big enough to make their experimental study realistic. The 0(a) 
order corrections to the lifetime of pionium are also calculated. 
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The pionium (7r + 7r~-dimesoatom) was firstly considered by Uretsky and Pal- 
frey |I|, although only after two decades an experimental study of such 
Coulomb dimesoatoms, initiated greatly by L. L. Nemenov |2j], became real. 
Recently about 200 pionium atoms, produced in a Ta target by 70 Gev pro- 
tons, were observed 0. But an accurate pionium lifetime measurement is 
still absent, and it is expected that this will be done in a more elaborate 
version of this j|] or other U proposed experiments. From this measure- 
ment a model independent information about pion scattering lengths can be 
extracted Q, 0,0. 

Already in [|T| it was realized that pionium can be formed in the elementary 
particle decays. Among them the K + — > 7r + ' A^ decay is most promising j3j. 
The branching ratios of this and other similar atomic decays are very small, 
making their experimental study an elaborate task. In the case of pionium, 
matters are further complicated by the fact that it is a very short lived object. 

Let us note that the only atomic decays of mesons, experimentally studied 
at yet, are 7r° — > 7^2e an d Kl - ► vA^ |jK)f| . Their measured branching 



ratios agrees well to what is expected theoretically [0, [TJ, |T2fl . Even 0(a) 
corrections for them were calculated [T3], [TJf, although some part of 0(a) 



corrections for Ki — ► vA^-, originated from the normalization constant, was 
omitted in Qnj . 



it is expected that several mesonic factories with very high luminosity 
will begin to operate in a near future. The study of atomic decays, even 
with pionium in a final state, can become real at these factories. Below we 
will consider some atomic decays of the elementary particles which can be 
interesting in this respect. 

Further information about dimesoatoms, not presented in this article, can 
be found in [jT5|, [16], [T7[ . 



2 Pionium in the 77, 77/ and K-meson decays 

In the nonrelativistic approximation pionium state vector coincides (up to 
irrelevent phase factor) to the following superposition jE|] : 



IP a \ f dqP A ; + ,,Pa - ^ U; m 



where ^(if ) is a momentum space wave function for A^ Coulomb atom and 



m = ■ s /m^+ « J Hf> appears because of relativistic normalization 



(?\t) = 2E p (2ir) d 6(?-?) , E p = \lT + mn? , (2) 

which we will use throughout this paper for any one-particle (composed or 
elementary) state vector. 

(1) means that an amplitude for the atomic decay P\ — ► P2 + A2 n is given 
by the following expression 

(PAIfi) = Jffli, % ^, Qi f^ 0) , (3) 



where ^(1? = 0) is a Schrodinger wave function for a hydrogenlike atom 
at the origin, and I(Qi,p + ,p-, Q2) is an amplitude for the nonatomic decay 
Pi(Qi) — > 7T+ (p+) + 71 " ~(p~) + ^2(^2)- If this amplitude is known, (3) allows 
us to calculate a ratio of atomic and nonatomic decay widths. 

If Pi and P2 are pseudoscalar mesons, then in a Pi-meson rest frame 
I(Qi,p + ,P-, Q2) depends only on two independent energies, and it is conve- 
nient to use instead of them the conventional Dalitz variables: 

V3(E + - EJ) M 1 -E + -E.-M l 

M 1 -2m-M 2 ' Mi -2m -M 2 ' y ' 



If M\ — 2m — M 2 <C Mi, it is expected [jTSj] that higher order terms in a 
power series expansion of I(x, y) will be suppressed and in a good approxi- 
mation 

\I(x, y) \ 2 ~ 1 + ay + by 2 + ex 2 

Note that terms linear in x are forbidden by CP-invariance. After integrating 
over two and three particle phase spaces: 

r(Pi -► p 2 ^) = 



/ KiV^lP!)! 2 -— (27r) 4 £(Qi - Q 2 - Pa) 



2Mi v y v ^ ^ y (2tt) 3 2P 2 (2tt) 3 2P a 

r(Pl -► P27T+7T-) = 

= / i/w., p+ , p -,q 2 ) i^moi - ^-^-^-) (2 rf | 8 S5_ 



3 m 3 



and using |\I/( x = 0)| = ^rp-, we obtain 



IXPi -> P 2 >W) 7T 3 TO 2 



= ««*(— n/(£, 



"-£ • (5) 



r(Pt -► P 2 7r+7r-) R v Mi 7 ' v ' w/l \ M{ 

x and y are Dalitz variables, which correspond to the atomic decay: 

X = ° V = 3 M 1 -2m-M 2 E2 = W, ■ (6) 

R is a dimensionalless remnant of the three-particle phase space integral: 

V ^+Jmax 7 fv^~- /max 7 | T / ^ |0 



R — dx + / dx_ |7(x + ,x_)| , (7) 

l*£+/min ^ (^ — Jmin 



£+ E 



where x + = j-^,X- = p-, and the integration limits in (7) are given by 



m_ 1 M 2 (2m + M 2 ; 



, 1 , 2m 2 -M| 

, X_ ) mm ; jTl ( 1 X -L 1(1 ZI_L H~ 7—7; ) 717 

^ ^m^ 2 (l-2x+ + f|) U A M? ; 



=F' 



x2 ^-)(l 2x 1 M 2 (2m-M 2 ) M 2 (2m + M 2 ) 



At last 



1, 



m 2 Ml 



(9) 



Ml Ml' ' 
The Dalitz-plot distribution for 77 — > 7t + 7t^7t° decay had been measured 



H| and the result is 



\I v (x, y)\ 2 ~ 1 - (1.08 ± 0.014)2/ + (0.03 ± 0.03)t/ 2 + (0.05 ± 0.03)x 2 (10) 

Inserting this into (5) and (7), we get (it is assumed that A^-k is produced in a 
IS state. If we sum up over all nS states, the result will increase E^Li 3«1.2 
times): 



T(r/ 



TTTT TT 



0' 



0.91 ■ 10~ 7 



;n) 



Let us note, that a value 3.9 • 10~ 7 , cited in |n|, corresponds to a theoretical 
prediction from effective chiral lagrangian |2U| I(x, y) ~ 1 — 0.55?/ and seems 
to be too optimistic, though the accuracy of quadratic terms determination 
in (10) allows, in principle, to increase (11) several times. 

The results, analogous to (10), exist for K + — » tt + tt + tt~ decay [|21j| : 



for K L 



\I K +(x, y)\ 2 ~ 1 + (0.2814 ± 0.0082)?/ - (0.001 ± 0.023)?/ 
-(0.099 ±0.019)x 2 , 

7T + 7r~7r decay 



Ifefe ?/)! 2 ~ 1 ~ (°- 917 ± 0.013)2/ + (0.149 ± 0.013)?/ 2 + 



+ (0.055 ±0.010)x 2 , 
and for 7/ — ► 77tt + tt~ decay [p3fl : 

l^/(^,2/)| 2 ~ |1- (0.08±0.03)2/)p 
Using them, we get 



r(if+ -> 7r+^ 27r ) 

T(i^+ -> TT+TT-TT ) 
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r(#i 



7r°A 27r ; 



I\i^ -> TT+TT-TT ) 



8.6 • 10~ 7 



r(r/ / -> 7/A 27r ) 
r(r/ 



1.4 -10 -6 (12) 



T/TT^TT J 

For i^ + — » TT+7427T decay, it is necessary to take into account the identity 
of Tr + -mesons, which increases the result two times. An extra ~ 5 times 



difference between K + and Kl decays is due to Ijc+(x,y)/lK L (x,y) ~ 2.15. 
Note that K + — > 7r + 74 27r decay was considered earlier in |8[ with slightly 
different result. 

Taking nonatomic decays branching ratios from [|23fl , (11) and (12) can be 
transformed to 

Br{r] -> 7r°A 27r ) « 2 • 10" 8 £r(7/ -> r/,4 27r ) « 6.2 • 10~ 7 

Br(K + -> 7r + A 27r ) « 5.5 • 10~ 7 £r(#£ -> 7r°A 27r ) « 1.1 • 10~ 7 (13) 

For 0- factory (13) means about 10 4 K-meson atomic decays with pionium 
per year. So, we think, the study of such atomic decays at 0-factory is not 
only realistic, but a desirable task. 

3 Pionium in the ip and T-meson decays 

For c - r and B-factories ij)(2S) -> ip(lS)A 2ir and T(2S) -* T(1S)A 2 „ de- 
cays can be interesting since the corresponding nonatomic decays have large 
branching ratios. 

In the nonrelativistic approximation the most general form for the V\ —> 
V 2 tt + 7t^ decay amplitude, which follows from PCAC, is [|25| 



(y 2 7r + 7r-|yi> = el • e 2 [-Aq + ■ q_ + B E + E_] + 

+C (el • q[ e 2 • q 2 + e[ ■ q 2 e 2 ■ q~i) (14) 

where A,B,C are approximately constant and ei,e 2 are Vi,V 2 vector meson 
polarization vectors. 

There are some theoretical indications [^5j and experiment confirms [)27 



that C — 0. Then to calculate the (5) ratio, we only need B/A ratio and it can 
be extracted from the tt + tt~ invariant mass distribution in the V\ — » V 2 7r + 7r~ 
decay. The results are |j27] : 

- = -0.21 ±0.01 (15) 

for the ip(2S) — > tl)(lS)A 27T decay, and 

B , s 

— = -0.154 ±0.019 (16) 

for the T(25) -> T{lS)A 2n one. 



After summing over vector meson polarizations, it follows from (14) that 
(if C=0) 

1 M 2 m 2 B 

(y 2 7r + 7r-|y 1 )--(l + ^|-2^2)-l + x + + x_--x + x_ , (17) 

x + and x_ were defined earlier. 

Using this instead of I(x,y) in (5) and (7), we get from (15) and (16) 

r(»/)(2S) -. j/>(lS)7T+7r-) ~ ' 
r(T(25) -> T(lg)^) _ 8 

r(T(2S) - T(lS)ir+ir-) 5 ° ' ( 8) 

which correspond to the following branching ratios 

Br(iJj{2S) -»- ^(15)A 2 ^) « 1.4 • 10~ 8 

Br(T(2S) -> T(lS)A 27r ) « 10" 8 (19) 

Unfortunately this is too small for B-factory. We can expect only several 
events per year. So it seems unrealistic to study T-meson atomic decays at 
B-factory. 

4 0(a) order corrections to the pionium lifetime 

The main decay mode for pionium is t4 27I - — > 7r°7r° and its amplitude according 
to Mandelstam J2?| can be expressed in the form 



ttVI 4*) = / -^- J(pi, ^ ^ + p, ^ - p) X (p) , (20) 



where x(p) is 74 27r -dimesoatom bound state Bethe-Salpeter wave function and 
J(pi,P2,p+,P-) stands for 7r + 7r~-irreducible kernel for the reaction 

7T + (p+) + 7T-(p_) -> TT ^) + 7T°(p 2 ) . 

Up to 0(a) terms, J is constant defined through pion scattering lengths 

ao and a 2 IM3 : 



32 

J = — 7rm(ao - 02) • (21) 

o 

So (20) can be rewritten as 

(A°\A 27T ) = J- X (x = 0) , (22) 



where 



dp 



x(fi = 0) = / T^uXip) 



(23) 



is a configuration space Bethe-Salpeter wave function at the origin. 

The Bethe-Salpeter equation for \(p) m a ^2n rest frame, up to 0(a) 
terms, takes the form 



2 \2 

m + p 



M A 



iX 



Tl 1 



dq 



2 

x(q) 



+ Po)' 



o \2 

m + p 



M A 



Po)' 



x(p) = 



(p ~ q) 2 ~ ^ e 



(24) 



where A = ^A 



(24) corresponds to the Wick-Cutkosky model |50f| . Let us note that this 
fact was firstly noticed and used to calculate 0(a) order corrections to the 
Kl —> vA^ decay width in [jll . 

for a review), a ground state (15 in the 



(25) 



According to [[3Dj (see also 
nonrelativistic limit) solution of (24) corresponds to 

-i g(z) dz 



x(p) = j_ x 



[A + Bz] 3 



where 



1 
4 



A = m 2 - -M\ -p 2 = A 2 -p 2 



B = p M A 



and g(z) spectral function satisfies the following integral equation 



9(z) = 



A a 



1 



If Ma = 2m — E, E <C m, then A 2 « mE and 

A 2 + \M\y* = W^E 6{V) 
because 



' Z S(z - y) + \±^B(y - z)]g(y) dy (26) 



(27) 



e 2 + y 2 
when e — >• 0. 

Substituting (27) into (26), we get 

Att 



7T(%) 



9\Z) = 



2mymE 



1 



1*1) 0(0) 



(28) 



So 

X * =1 
2mvmi? 

which really gives a hydrogenlike atom nonrelativistic ground state energy 
level E = m g-. 



But the solution (28), found in [|30j| , is not a complete 0(a) order solution. 



Indeed, taking g(z) = go(z) + agi(z) and 

1 2tt 



A 2 + \M\y 2 m 2 a 



6(y)+a(y) , (29) 



where a(y) has 0(a) order smallness compared to the first ~ 5(y) term, we 
get from (26) (Nq and Ni are constants) 

g (z) = N (l - \z\) 9l (z) = iVi(l - \z\) + A /' <jfo)#(*, i/)pu(y) dy (30) 
where 

#(*, y) = ^e(z - y ) + \^e(y - z) . 

\-y l+y 

Calculating in the a — > limit the integral in (30), we get 

9l (z) = N x (l - |*|) + ^{(1 - |*|) ln(«) + (1 + M)[m(2|*|) - ln(l + |*|)]} . 

Therefore, the complete 0(a) order solution of (26) looks like 

g(z) = N{(1 - \z\) + -(1 + |*|)[ln(2|*|) - ln(l + \z\)}} (31) 

Substituting this in (25), we get the 0(a) order pionium Bethe-Salpeter wave 
function 

X(P, p a) = 

= (A 2 -p 2 ) [m 2 - (f + p) 2 } [m 2 - (f -p) 2 } V + k Xi{p > Pa) ) (32) 



where 



Xi(P, ^) = m2 2(A ^ ^f ln(m 2 - (£* - pf) + 

m 2 - (f + 2?) 2 ln(m 2 _ (^ + p )2) _ i n (A 2 - p 2 ) + 0(« ln(a)) (33) 



2(A 2 -p 2 ) v v 2 



A normalization constant N is defined from the normalization condition 
32|1 , which in the 0(a) order takes the form (^tt rest frame is assumed) 

dM 2 



9 9 ' 

—2m a « a- 

zA 2 /• d£> 



But the integral in the r.h.s. of (34) equals to 



~2~ ~ Po) ~ m 2 ~ P 2 



(34) 



m d 2 r r dpp 1 1 

2d(A 2 ) 2 J P J 27riA 2 +p 2 -p 2 [ {^+ Po y-m 2 -p 2 
1 



(^ - Po) 2 -m 2 -f 

Remembering that in fact m 2 in the above expression should be replaced by 
m 2 — ie, we can perform an integration over dpo and obtain 

2l7T 2 roo x 2 1 4(m 2 + x 2 ) 3 

" MfX dx (A 2 + x 2 ) 2 Vm 2 + x 2( A 2 + x 2 )_ VA 2 + x 2 ^~ 



m 4 a 3 7r 



1-2- + OM) 



So (34) takes the form 

i 2 ™ 2 = 

(2tt) 4 m 4 a 3V tt 



^99 i ,.9 ^ 3 / 0K- 

2m 2 a 2 = - — -iV 2 — -1 - 2- 



and therefore 



A = 32 v / ^(-ma) 5/2 (l + -) . (35) 

2 7T 



Analogously 



dp i A roo x 2 2 2 

X\P) ~ 0_9 7\#2 n ^ -—-—jr{—j===== -==—= + 



(2tt) 4AUV 87r 2 Mi^o A 2 + x 2 Vm 2 + x 2 VA 2 + x 2 






^m 2 + x 2 )(A 2 + x 2 ) 167rm 2 a tt 

Therefore 

x(* = 0) = - ^ . (1 + -) » -i=*(* = 0) (1 + 2^) (36) 



Substituting this in (22), we finally get 



T(A 



2tt 



7T 7T , 



0„0n 



r (^^7r u 7r u )(l + 4 



a. 



7T 



(37) 



where 



0„0n 



T (A 27r -► 7r u 7r 



16tt 



ao — a<i, 



\ 



2(111^+ — m^ 



m T 



|#(£=0)f 



~ Xi part of the Bethe-Salpeter wave function contributes neither in its 
value at the origin nor in its normalization, because in the 0(a) order 



h 



adpt 



(^+P0 



m 2 — p 2 



M A 



PO 



m 2 — p 2 



i7ra5(po) dpo 
M A (A 2 +p 2 ) 



and 8(po)xi(p; Pa) = 0. 

As a last remark, let us note that a nonrelativistic approximation (3) for 
the Mi — » M.2 + A^ decay amplitude follows from 



and 



X{x = 0) 



m 



*(£ =0) 



P2 ^ |Pl> "/(2^ 
^(Qi,y,y,Q 2 Mx = 0) 



AQi> -y + P> y - P> Qt)x(p) 
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